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ABSTRACT 

A star on a nearly radial trajectory approaching a massive black hole (MBH) gets tidally disrupted 
if it comes sufficiently close to the MBH. Here we explore what happens to binary stars whose centers 
of mass approach the MBH on nearly radial orbits. The interaction with the MBH often leads to 
both stars being disrupted in sequence. We argue that such events could produce light curves that 
are substantially different from those of the single disruptions, with possible features such as two local 
maxima. Tidal forces from the MBH can also lead the binary components to collide; these merger 
products can form highly magnetized stars, whose subsequent tidal disruption may enable prompt jet 
formation. 

Subject headings: binaries: close — galaxies: kinematics and dynamics — galaxies: nuclei 


1. INTRODUCTION 


Tidal disruption of stars by massive black holes 
(MBHs) in galactic nuclei has been studied by theorists 
for over 40 years, in the co ntext of supplying gas for ac¬ 
cretion on to SMBHs (e.g., |Hillslll975 |Rees |1988| jPhin- 
ney||1989D. During the past sixteen years, several tidal 
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Observational 


progress has stimulated ongoing theoretical efforts to un¬ 
derstand the details of tidal disrup tion physics and to 
model the lightcurves of TDEs (e.g., Lodato et al.||2009j 
Strubbe fc Quataert||2009[ |G ui llochon fe Ramirez-Ruiz| 
20131120151 IShen fc MatznCT|20i4f]5Hiokawa et al.||2015[ 


Bonnerot et al. 2015p. The advent of deeper all-sky sur- 

veys such as LSS'i', as well as forthcoming transient sur¬ 
veys in the radio and ultra-violet bands, are anticipated 
to considerably increase the number of TDE observa¬ 
tions. 

The majority of the stars in the field, however, are in 
binaries, and some of the stars on a collisional course with 
the MBH will be members of a binary. It is therefore 
interesting to consider an encounter between an MBH 
and a stellar binary, whose center of mass is approaching 
the MBH on a nearly parabolic orbit with the pericenter 
distance on the order of the tidal disruption radius of 
a single star. In this paper we show that such encoun¬ 
ters often lead to the tidal disruption of both stars in 
sequence. We find that when the encounter between a 
stellar binary and an MBH leads to a TDE, both mem¬ 
bers of the binary are disrupted in sequence more than 
75% of the time. We estimate the frequency of such dou¬ 
ble TDEs, and argue that many of them will produce 
light curves that are different from those of single TDEs. 
We discuss the implications of this scenario for recent ob- 
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servations of double flares in SDSS J015957.64-1-003310.5 
( Merloni et al.|[2M5 ) and IC 3599 (Grupe et al 2015). 
Collisions betweeii binary componen ts caused by tidal 


forces from the MBH are als o common (Ginsburg & Loeb 
20071 jAntonini et al.||2010|). They populate the phase 


space of low angular momentum stars with merger prod¬ 
ucts many of which will eventually be tidally disrupted. 
These merger products will have much larger magnetic 
fields than single stars (Wickramasinghe et al.|2014 Zhu 
et al.|2015 ). The high magnetic flux is accreted onto the 
iVIB H during a TDE and can aid in generating prompt 
jets (Tchekhovskoy et al. 2014), such as observed in Swift 
J164449.3+573451. -- 

The structure of the paper is as follows. In sectionj^we 
develop a population model for the binaries approaching 
the MBH with small impact parameters, analytically jus¬ 
tify why there should be a significant fraction of double 
TDEs, and describe the setup of numerical experiments. 
In section]^ we present numerical experiments in which 
we follow the Newtonian dynamics of stars in binaries 
approaching the black hole, and flag tidal disruptions of 
individual stars and stellar collisions. In section |4] we 
discuss the observational implications of our findings. 

2. MODEL OF THE BINARY POPULATION 

We are interested in binaries whose nearly radial ap¬ 
proach to an MBH leads to a possible disruption of one or 
both binary companions. These binaries become tidally 
separated before they reach the point of closest approach 
( [HUIs|[TM| ISari et al.||2010| ). 

A star 01 mass m* and radius i?* will be tidally dis¬ 
rupted by an MBH of mass Mbh if its orbit brings it 
within the tidal disruption radius of the MBH, 


7?*td = 


( 1 ) 


here Q ~ 1 depends on the structure of th e star ( Phinney] 


1989 


Guillochon & Ramirez-Ruiz 2013). From here on 
all assume a = 1; however, some sur viving stars 


we s. 

will, in fact, lose a fr action of their mass ([Guillochon 


& Ramirez-Ruiz 2013). A binary on a nearly parabolic 


orbit around the iVIBH whose center of mass has a peri- 
apsis Tp < i?*TD is likely to first be tidally separated, and 
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Fig. 1. — An illustration of loss-cone dynamics. The red rectan¬ 
gular region shows the angular momenta of individual stars which 
will be tidally disrupted, provided nothing changes their orbit on a 
dynamical timescale. The blue dashed lines indicate the loss cone 
for tidal separation of stellar binaries (not to scale). At small radial 
separations, r < Tmim the binary loss cone is empty, so binaries 
like the one on the bottom are tidally disrupted before entering 
the regime where individual stars can be disrupted; therefore, they 
cannot be the source of double TDEs. At larger radii r > Tmin 
(top binary), the loss cone is full and direct penetration into the 
stellar tidal disruption region is possible. Double TDEs originate 
from the full binary loss cone region. 

then to have both stars undergo sequential tidal disrup¬ 
tions. This happens because the component stars have 
specific angular momenta that are close to that of the 
binary’s center of mass, and thus undergo a similar pe- 
riapsis approach. The specific center of mass angular 
momentum of the binary on an orbit around the MBH 
with Tp = i?*TD is 


l*TD Y GMbkR* ^ • (2) 


The fractional difference in the angular momenta of the 
companions when the binary with semimajor axis a is 
tidally separated is of order 



Therefore, to produce a double TDE, a binary must en¬ 
ter the “red” region by jumping in from outside the blue 
lines! 

Consider this dynamics in more detail. Binaries wan¬ 
der stochastically in angular momentum space by gravi¬ 
tational scattering on other stars; this process can be ac¬ 
celerated by gravitational interaction with massive per- 
turbers outside of the MBH’s sphere of influence (SOI) 
( Perets et al.||2007 ) and by secular gravitational torques 
in the near-KepIerian po tential deep inside the SOI 
(Rauch & Tremaine 1996). If the typical change in an- 
gular momentum after one orbit is Al <C ?lCi the binary 
can never get onto a small angular momentum orbit in¬ 
tact. In the illustration of Fig. the lower binary will 
be tidally separated immediately upon crossing the dot¬ 
ted blue line and will never reach the red r ectangle; this 
regime is kno wn as the empty loss cone (Lightman &:|| 
Shapiro 1977). Meanwhile, if Al > ZlCj the binary can 
take large steps in angular momentum space, and can im¬ 
mediately jump anywhere within the loss cone I < Zlc! 
this is more probable for binaries coming from far aww 
from the MBH, as illustrated by the top binary of Fig.[lj 
In the full loss cone regime, any angular momentum 
within the loss cone is nearly equally likely, and the prob¬ 
ability that the binary ends up with a two-dimensional 
angular momentum < |/*td| is ~ (^*td/^lc)^- 

The dividing line between the full and empty loss cone 
corresponds to Al ^ /lc ■ If the relaxation is governed by 
two-body interaction with other stars, the timescale for 
the angular momentum to change by the specific angular 
momentum of a circular orbit, Zcirc = a/GITbht, is 


0.065U'’' 


log A ’ 


( 5 ) 


where n{r) the local stellar number density and v ^ 
^/G' MRH7r the typical velo city at distance r from the 
MBH ( Spitzer fc IIart|[l97I ). Assuming the stellar num¬ 
ber density is a power law n(r) = no(r/ro)'’', the loss 


cone is full for r > Tmin with 


7/3 


AfsH ) 
m* J 27rnologA" 


O.OGSr)) 


1/(4+7) 


( 6 ) 


At r > Tmin, the ionization of bi naries through interac- 
tions with other stars is unlikely (Hopman 2009). 

The rate at which binaries enter the region where in¬ 
dividual tidal disruptions are possible is roughly 


i.e., less than one for m* = Mq, R* = i?©, Mrh = 
10®Mq and a < 50 AU. Therefore, the companion peri- 
apses will also typically be similar. 

The phase space where angular momenta are suffi¬ 
ciently low to allow for stellar tidal disruptions is marked 
by a solid red rectangle in Fig. It is non-trivial for a 
binary to be dynamically inserts into this region. This 
is because binaries will be tidally separated by the MBH 
whenever the binary’s angular momentum around the 
MBH is less than 


d^N 

drdadt 


(a,r) 



^ ^ ^min 1 

otherwise. 


where we take the density function of binar y semi-maio r 
axes to be 7r(a) = [ln(amoa:/amm)]“^(l/a) ( |6pik||l924[ ), 
and 


P{r) = 27r 



( 8 ) 
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a region bounded by the dotted blue lines in Fig. 


is the period of the center of mass of the binary around 
the MBH. 

The loss cone is “half-full, half-empty”: roughly com¬ 
parable numbers of binaries should enter it in the full and 
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loss-cone dynamics in galactic nuclei). 
Therefore, the rate of mergers per unit semimajor axis 
per unit time, i.e., the integral of Eq. over r, can 
be roughly approximated by evaluating the integrand at 
T = Train and multiplying by This yields the over¬ 

all distribution of semimajor axes among binaries whose 
components can undergo individual TDEs immediately 
following the tidal separation of the binary: 


(PN 

dadt 


0, a-7/(8+27) ^ „-14/9 ^ 


(9) 


where we assumed a Bahcall-Wolf cusp, 7 = 
Shal lower cusps, such a s the one in the Galactic cen¬ 
ter (|Bartko et al. 2010), would increase the fraction of 
double I DEs among all TDEs. This estimate is simplis¬ 
tic, and does not take into account the detailed stellar 
distribution outside the SOI. However, a more soph is- 
ticated analysis performed by MacLeod et al. (2012) in 
order to study the rate of tidal disruptions of giant stars 
is in agreement with our simple scaling. Eor a broad class 
of models, they found that the probability per unit time 
for a giant to be disrupted scales as p{R*) oc i^o. 35 ±o.i 5 ^ 
Applying our calculation for the tidal separation of bi¬ 
naries, Eq. (p), to tidal disruptions, we would predict 

p(R* ) oc i?* , in good agreement with (MacLeod et al. 

20 T^. 

Further integrating Eq. ([^ over all a, we can obtain the 
total rate of binaries entering the stellar disruption loss 
cone and compare this to the rate of individual TDEs 
of single stars, which is estimated by integrating the 
same equation with a replaced by i?*. Assuming the 
same density n(r) of binary and single stars, and us¬ 
ing the distribution choices discussed below, as many 
as (9/5)[ln(amaa;/omm)]“Hamm/-R*)“®^® ^ 10 % of all 
TDEs could be contributed by double disruptions. 

We now describe the details of our simulated popu¬ 
lation. We fix the MBH mass to Mbh = 10®Mq in 
all cases. The primary compone nt mass is dra wn from 
the Kroupa initial mass function ( Kroupa|2001 ) and the 
secondary is drawn from a distribution of the mass ra¬ 
tio q = TO 2 /TO 1 given by p{q) oc for q G [ 0 . 2 , 1 ], 

broadly consistent with the observed field di stributions 
of bot h low-mass dPuguenno y fc Mayor||1991|) and high- 
mass (Sana et al.|2U12 ’21)13) stars. We cut off the mass 
distribution for q < U.2 and mi 2 < 0.1 Mq, where it is 
poorly observationally constrained. 

As discussed above, binaries of interest are distributed 
according to p(a) cx so we simulate p(Rb) cx 

Pfj for binary periods Pb bet ween 0.1 day and 1000 
years. Approximately following [Duquennoy fc Mayor] 
(1991), we assume that binaries with periods less than 10 
days are tidally circularized, binaries with periods above 
1000 days have eccentricities drawn from a thermal dis¬ 
tribution p{e) = 2 e, and binaries with periods between 
10 and 1000 days have their eccentricity drawn from a 
Gaussian distribution with mean 0.3 and standard devi¬ 
ation 0.15, cut off at 0 and 1. We impose the additional 
constraint that the initial binary periapsis must be at 
least twice the sum of their radii, where the stellar radii 
are i?* = with k = 0.8 for m* < M© and 


k = 0.6 for m* > M© (|Kippenhahn & Weigert||1994 ). 

We start each binary on a parabolic center-ot-mass tra¬ 
jectory around the MBH at a distance ten times the bi¬ 
nary tidal separation radius, with MBH periapses drawn 
uniformly from G [0,300] i?© (Fig. demonstrates that 
no TDEs happen for larger periapses, in line with ex¬ 
pectations, Eq. 0 )- The three-b ody evolution is carried 
out with the RE BOUND integrator (Rein & Liu 2012 Rein 
fc Spiegel||2015[ ). The evolution proceeds until either (i) 
the components collide (the distance between the stars 
becomes smaller than the sum of their radii; we label all 
such events as ‘mergers’ be low although individ ual stars 
may survive in some cases ( [Antonini et al.|[2011 )); or (ii) 
one of the components is tidally disrupted (comes within 
a distance i?*TD of the MBH), after which the future 
trajectory of the remaining companion can be computed 
analytically as a 2 -body problem; or (hi) the binary com¬ 
panions fly past the MBH without being disrupted or 
merging. 

3. RESULTS OF NUMERICAL EXPERIMENTS 


We numerically simulated the dynamics of 1000 bina¬ 
ries as described above. Among these, 18% resulted in 
sequential tidal disruptions of both binary components 
and 5% resulted in only one component (typically the 
more massive one) being disrupted. Additionally, binary 
components merged in 6 % of simulations before either 
was disrupted; more mergers will arise from larger initial 
periapsis radii than considered here. The remaining 71% 
of simulations produced no disruptions or mergers. 

Figure shows the relative frequency of the possible 
outcomes as a function of the center-of-mass periapsis 
around the MBH and the initial binary semi-major axis. 
Both components are likely to be disrupted if the peri¬ 
apsis is sufficiently small; see discussion around Eq. (|^. 
Single TDEs are relatively more common in comparison 
to double TDEs for large periapses, where they are corre¬ 
lated with significant companion mass differences which 
lead to different tidal disruption radii. Mergers are rela¬ 
tively common for small binary separations, as tidal per¬ 
turbations to the orbits are more likely to lead to col¬ 
lisions between companions; mergers preferentially hap¬ 
pen for larger periapses, where disruptions are avoided. 

When both components are tidally disrupted, we can 
ask whether individual events are resolvable. The differ¬ 
ence in peak luminosity times, ATpeak, is the difference 
in the rise times tpeak to the maximum value of the ac¬ 


cr etion rate Mr,Bak, estimated using t he fitting formulae 
of Guillochon & Ramirez-Ruiz (2013), corrected for the 
the ditterence in the times when the binary components 
reach their periapses. An additional correction of order 
Rsf /a to peak times and mass accretion rates, due to the 
energy shift of individual stellar orbits away from the 
zero-energy parabolic orbit, is not included. Figure 
shows a scatter plot of the ratio of ATpeak to the smaller 
of the two rise times, against the ratio of peak accretion 
rates for the two disrupted companions. We see that 
peak accretion rates are typically comparable (the ratios 
are within a factor of two in > 70% of cases), and never 
exceed 10 : 1. For 20% of double TDEs, peak accretion 
times are separated by at least half the smaller rise time; 
even larger time delays ATpeak > 2 months with indi¬ 
vidual rise times of tpeak < 1 month are observed. The 
eccentricity vectors of the two component orbits, which 
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Fig. 2.— Distribution of mergers, double disruptions, single 
disruptions, and uneventful fly-bys (from top down) among the 
1000 numerical simulations drawn from the distribution described 
in Sectionl^ as a function of the center-of-mass periapsis Vp around 
the MBH (top) and binary semi-major axis a (bottom). 


point to the periapses, are nearly perfectly aligned in all 
cases, but orbital plane misalignments are occasionally 
possible, though rare (only 1 of 183 simulated double 
TDEs has orbital planes misaligned by more than 90 de¬ 
grees) and generally correlated with larger AT. 

As an illustrative example, we consider a binary with 
O. 4 OM 0 and 0.27Mq components from our simulations. 
The initial semimajor axis is S.2Rq and the initial center- 
of-mass periapsis radius is 42i?Q. The primary is dis¬ 
rupted first, but given the small binary separation, the 
secondary follows on nearly the same trajectory and is 
disrupted a few minutes later. The secondary has a peak 
rise time of 46 days, compared to only 23 days for the 
primary. We show an illustration of this double TDE in 
Fig.R 

When there is a single TDE, the surviving companion 
remains on a small angular momentum orbit. Half of 
the surviving companions are on bound orbits, with a 
median return period of ~ 50 years. The surviving com¬ 
panion may undergo a partial tidal disruption on the first 
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Fig. 3.— Double TDEs: the ratio of the difference between 
peak times and the shorter rise time, vs. the ratio of the peak 
mass accretion rates. The color (red is 80 days, blue is 20 days) 
represents the shorter rise time. 



Fig. 4.— An illustration of the mass accretion rate for a typ¬ 
ical double TDE described in the text, const ructed using fitting 
formulae of|Guillochon &: Ramirez-Ruiz| ||2013[l. 


periapsis passage, and be subsequently disrupted after 
returning to the MBH. 

Binary collisions fall into two classes. Some 30% of 
mergers happen en route to the periapsis; the merger 
product arrives at periapsis after no more than two nom¬ 
inal dynamical timescales of the merged star, so we ex¬ 
pect such merger products to be tidally disrupted even 
if the MBH periapsis is a few times the nom inal tidal 
disruption radius of the merged star (see also |Antonini| 
et al.|20fT). The other 70% of collisions happen after the 


stars pass periapsis, sometimes tens of tidal disruption 
radii away from the MBH; the close approach of binary 
co mponents after per iapsis passage had been pointed out 
by Sari et al. (2010). These merger products have low 
angular momenta around the MBH, and most are likely 
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to be tidally disrupted on future orbits as they undergo 
a random walk in angular momentum space. 


4. DISCUSSION 

We have carried out dynamics simulations of tidal in¬ 
teractions of stellar binaries with an MBH. We found 
that both binary components are tidally disrupted in 
18% of our simulations. Double TDE events may rep¬ 
resent nearly 10 % of all stellar tidal disruptions. 

Double TDEs may have a distinctive d ouble-peak sig¬ 
nature. Intriguingly, Merloni et al. (20151 recently inter¬ 
preted the flare in AGIN bDbb 4015957.64-1-003310.5 as 
a possible TDE, but with an unexplained second peak 
in the light curve. A prompt double TDE could provide 
a similar signal. However, while Eig. can be viewed 
as an illustration of a double-peak light curve, we stress 
that we have not simulated the hydrodynamics of such 
an event. Such simulations, including the possibility of 
orbital misalignments between tidal streams which we 
observe in some cases, are necessary to predict the full 
signature of double TDEs. 

A smaller fraction of our simulations resulted in single 
TDEs. In half of these cases, the surviving star is left 
on a bound orbit with a period ranging from 6 months 
upwards, and a median period of half a century. If the 
surviving star was partially tidally stripped during the 
initial encounter, it could lose more mass on a subsequent 
periapsis passage and produce a secondary flare. The 
recent observation in IC 3599 of a repeated flare 20 years 


after a ROS AT event initially classified as a TDE (Grupe 


et al. 


20151 could be representative of this scenario. 


Stellar collisions after periapsis passage can leave 
merger products on low angular momentum orbits 
around the MBH. These merg er products are likely 


to ha ve strong magnetic fields (Wickramasinghe et al. 


20141. Such magnetic fields could play a critica l role in 


powering prompt jets (Giannios fc Metzger 2011|) such as 


observed in Swift J164449.3-1-573451 ( Bloom et al.||2011 


Burrows et al. ]|20111 ILevan et ap|20ll[ [Zauderer et al. 


2011[ ). Tchekhovskoy et al. ( |2014p argue that a magnetic 
flux > 10^^ G cm^ is required for prompt jet formation 
but is unlikely to be delivered by the star itself, and must 
be generated by the accretion-disc dynamo or be already 
present in the MBH vicinity. However, the magnetic field 
in the merger product may ho st 0 . 2 % of the remnant’s 
total energy ( Zhu et aL]|2015 ), which implies a flux of 
^ 0.07 -\/Gto* = 4 X lO^®(m*/M 0 ) G cm^. This may 
serve as a useful seed for generating the field required for 
launching a jet. 
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